How fast can fluids squeeze through micro-pores? 
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We use a one dimensional symmetric exclusion model to study pressure and osmosis driven flows 
through molecular-sized channels, such as biological membrane channels and zeolite pores. Ana- 
lytic expressions are found for the steady-state flow which show rich nonlinear behavior. We find a 
maximum in the flow rates depending upon pore radius, pore energetics, reservoir temperature, and 
driving force. We also present exact mean-field results of transport through pores with internal de- 
fects. The interesting nonlinear dependences suggest numerous diagnostic experiments in biological 
and zeolitic systems which may reveal the features presented. 



The answer to the title question is tremendously im- 
portant in analyzing biological and industrial processes, 
and has received recent attention with the finding that 
tracer particle motions in a micro-pore are governed by 
subdifFusive dynamics Biological examples include 
membranes containing molecular-sized channels specific 
to water transport which participate in cellular volume 
rcgrilation controlled by hydrostatic or osmotic pressure 
[^,||. Classes of biological ion channels have also been 
demonstrated to be pores of molecular size Man- 
made materials such as zeolites and carbon nanotubes 
also contain many microscopic, statistically nearly single- 
file channels which can selectively absorb fiuids. This size 
specificity can be exploited in separation of a mixture of 
linear and branched chain alkanes, where the zeolite acts 
as like a sponge absorbing only the desired specie(s) [Q. 
Confining particles in zeolite pores can also serve to cat- 
alyze reactions: How fast can one get reagents into micro- 
pores and the products out? Therefore, the design and 
manufacture of porous materials is an economically mo- 
tivated area of research. Not surprisingly, vast amounts 
of numerical simulation have been performed on a variety 
of specific systems [Q. 

Confined particles are strongly interacting due to ex- 
cluded volume as shown in recent NMR and theoreti- 
cal studies Anomalies in numerically computed (MD) 
"diffusion" constants have been found However, nu- 
merical simulations neither access the long time scales 
required to study steady state fiow, nor offer a unifying 
physical picture of the parameter regimes important for 
transport. To obtain reasonable fiow rates using MD, ar- 
tificial external forces such as gravity are often imposed 
in ID systems such external forces can yield qual- 
itatively different behavior (such as shock profiles) from 
osmosis and pressure driven fiow , which occur in the 
absence of such intrinsic forces. 

A model of micro-flow that physically describes mi- 
croscopic transport and how it depends on macroscopic 
experimental parameters would serve as a useful bench- 
mark in more sophisticated models and complement more 
detailed MD simulations. Consider the molecular-sized 
pore shown in Fig. |l], with particles driven from [L) 
to (i?) either by osmotic "pressure" AH, or by hydro- 



static pressure AP. The pore is divided into sections i 
of length £. Enthalpic energy differences between bulk 
and sectioned pore particles are shown, as well as pos- 
sible activation barriers. Entrance(exit) rates at the left 
and right boundary sites are denoted by 0(7) and 5{I3) 
respectively. Here, adt and 5dt are the probabilities for 
pore entry in time dt only if the occupations Ti=i = 
and Ti^N = respectively. The probability per unit time 
a randomly picked interior particle moves one section to 
the right(left) is p{q) only if the site to the right(left) is 
unoccupied. 
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FIG. 1. Schematic of osmosis and pressure driven flow 
through membrane pores separating inflnite reservoirs (L) 
and (R). The coefficients a,f3,"f,5 are conditional solvent 
entrance and exit probabilities at pore ends. 

The assumption of no pass pores is accurate even for 
pores with diameters few times that of the particles: 
Overtaking requires a restricted subset of geometries and 
will be statistically rare. Even when overtaking occurs, 
interchange between particles at sites i and i + l does not 
contribute to net flux. Although our treatment is rigor- 
ous for no pass pores, it is accurate as long as the average 
number of particles in a cylindrical slice of length ^ is ^ 1. 
Consider the instantaneous number flux between sections 
i and i -|- 1, 
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J{t) = pT,{t){i - n+i{t)) - qn+i{t){i - n{t)) 
= p(n(t) - n+i(t)) + en+i(t) ~ eTi(t)Ti+i(t). 



(1) 



Although much attention has focussed on the asymmet- 
ric exclusion model, particularly in the thermodynamic 
limit, where many exact results are known ||l0|| , the sym- 
metric exclusion model is valid in the absence of external 
electric or gravitational forces, e = p — q = Q. This re- 
veals that hydraulic and osmotic pressure driven flows 
through microscopic pores are simply a consequence of a 
density gradient which induces particles to directionally 
or cooperatively diffuse; however, locally, a particle is as 
likely to move to the left or right if both left and right 
adjacent sites are unoccupied. The random updating in 
particle motions implicit in the model dynamics is clearly 
valid in systems where long wavelength collective modes 
are irrelevant. 

Linearity of J{t) when e = renders the mean-field 
steady state current J found from time averaging Eq. 
(|l|) exact. Steady state particle conservation along the 
pore length results in a linear density profile J — p{ti — 
tn)/{N — 1); this, along with the steady state boundary 
conditions Ji = J — Q!(1 — ti)— 7x1 and J^r — J = /Jtat — 
i5(l— r^v), determine the steady state particle number flux 
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p{a[3 — 7(5) 



(A^- l)(a + 7)(/3 + 5) +p(a + /3 + 7 + (5) 
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The kinetic parameters {/z} = (a, /?, 7, 6) are related to 
the relative enthalpies of activation between pore and 
bath particles. We assume local thermodynamic equi- 
librium (LTE), particularly valid in liquid phase osmo- 
sis across single biological pores where J < 10^/s, typ- 
ical pore diameters and interparticle spacings A ~ 5A, 
and ambient thermal velocities — 4 x lO^cm/s, yield 
collision times tcoii — X/vt — Ips <^ J^^. Therefore, 
particles suffer O(10'^) collisions before they are osmoti- 
cally transported, sufhcient for (LTE). As an illustrative 
example, we consider an axisymmetric right cylindrical 
pore. The kinetic parameters under LTE will be defined 
by simple transport theory, 

p « {vt/I) expi-Ep/ksT)- f3 w {vt/£) expi-Ep/keT) 

a = aoexp^-Ea/ksT) « InLVriTrrp) exp{~Ea/kBT) 

(3) 



where vt ~ y/k^Tjrn is the thermal velocity, I is cho- 
sen to be the minimum statistical spacing between pore 
particles, and is the number density in the left reser- 
voir. The minimum particle spacing i can be estimated 
for each set of {/i} from a ID Tonk's gas Since we 
explore considerable variations in {/i}, we choose £ ^ a, 
approximately a hard core repulsive diameter; thus, p 
given in (j^) represents a ballistic travel time over the 
distance £ ~ a, weighted by an energetic binding Ep. 



A larger choice for £ can be made if {£/a)Ti < 1 and 
with p appropriately rescaled and entropic factors in- 
cluded for {fi} (the E^ are then effective free energies); 
this is useful in multiple species models where steady 
state flows for long pores cannot be obtained analytically 
]l2[ . Note that for ^ a, local free diffusive transport 
described by p ~ £~'^ may obtain. But when € ~ a, 
E^ represent enthalpies determined entirely by molecu- 
lar potentials. Thus, Ea{r < Vp) — Ea{Q) ^ fcsT (where 
Ea(r = 0) = Ea) defines an effective pore radius r^. For 
pores that repel particles (top curve in Fig. |l|) and have 
negligible activation energies Ep and Ep, p/(3 ^ 0(1). 

Further simplification is gained by considering micro- 
scopically symmetric pores where /3 = 7, and normaliz- 
ing (denoted by an overbar) all quantities by P (such that 
J = 1 would be the maximum flow rate possible, when 
Tat = 1) so that 
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where 



apA 



{N - l){a + l){a + l-aA)+ p{2a + 2-aA)' 

(4) 
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dute ^AE/kBT 
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and AE = Ea - Es{Pr -Pl)- Eqn. (1) represent differ- 
ences in number density and/or enthalpies between (i?) 
and (L) and along with (^) determine the flow through 
symmetric pores. In osmotic flow under isobaric condi- 
tions, AE = 0, in pressure driven flow of nearly ideal 
gases, {dAE/dPii)T — 0, while pressure driven flows of 
liquids is described by [dnR/dP)T — 0. In the first two 
cases, the flux is predominately the consequence of in- 
creased permeable particle density in one of the reser- 
voirs over the other, while pressure driven flow of liquids 
results mainly from the relative reduction of pore en- 
trance activation energies brought about by hydrostatic 
compression. We first consider fixed A = 0.02 which 
corresponds to an osmotic pressure in aqueous solution 
of An ~ 25 atm or a hydrostatic pressure difference of 
AP ~ 0.025 atm of gas at STP. Using the Maxwell rela- 
tionship for molar volume, —{dAE/dPiijT — v, we find 
A — 0.02 also corresponds to Pr — ~ 25atm in pres- 
sure driven flow of water at 300K. 

When aA/{a + 1) is negligible, flows are essentially 
linear in A and defined by hydraulic or osmotic perme- 
abilities, J = LpAP or J — Pos An. Although Lp and 
Pos have often been, and continue to be interpreted using 
macroscopic fluid mechanics a microscopic descrip- 
tion arises here. In the limit p ^ {a + 1)N, rate limiting 
steps involve pore entrance or exit. Linearizing (|^) and 
(0), we find 



a / aA^ \ 
2{a + l)kBT \ dPR J , 



(6) 
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for pressure driven flow of dense liquids. The tem- 
perature dependence of Lp will be determined by 
-Efs/ksT [-Ec^/ksT] for a > 1 [a < 1] if the thermal 
coefficient of expansion « 0. When p ^ {a + 1)N, 



aq 



f d/\E \ 



(7) 



which has a {Ec,-Ep~Ep)/kBT [{Ep ~ E^ - Ep)/kBT] 
Arrhenius temperature dependence for a ^ 1 [a ^ 1]. 
Note that the a ^ 1 regime of (0) yields a curious 
Lp oc dependence. In the limit where holds, 
the bottlenecks occur in the particle motions within the 
pore interior. Expressions for Lp of ideal gases and 
Pos{Pl — Pr) are found from ^ and (0) by replacing 
(dAE / dPii)T respectively in the corresponding 

limits; the temperature dependences remain unchanged. 
Values of all energetic barriers mentioned above are qual- 
itatively consistent with osmotic flow through biological 
pores where hydrogen bonds {Ea ~ lOfc^T in the a ^ 1 
limit of (0)) must be broken before water can enter. 
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FIG. 2. (a) 1.5 X 10^ J{N, a) for A = 0.02 and fixed 
po = 1.0. Various lengths are indicated. The solid curve 
between and 1 is the average occupation number of any 
site. On this scale, the diflference n — rjv is not apparent but 
varies qualitatively as J. (b) The solid curve is 10'' X J(5) for 
fixed Po = 0.1 See text for explanation of curves P0P1P2 and 
PoPiPiPs- 

The solid curve in Figure ||(a) shows (1.5 x 10^) x J{N) 
(from Eq. (|^)) for p = 1, A = 0.02 for various pore 
lengths N ~ L/a. We find a maximum flux J* at 



2po + (V - 1) 



(N-m-A) 



1/2 



(8) 



for fixed p — pQ. The maximum occurs because at in- 
termediate occupancies, the pore is conducting a sig- 
nificant number of particles, but is not yet choked off. 



Large values of a represent pores which are attractive 
to the solvent, depicted by the dashed energy landscape 
in Fig. ^ As the pore is made increasingly attrac- 
tive Ep > 0, and /3 must eventually diminish so that 
p — p/f3 oc e'>q)[{E/3 — Ep)/kBT]. However, we will 
explicitly show that the maximum at J* can persist. 
Assume no activation barriers at the pore mouths, ie., 
Ep{Ea) — for EaiEp) > and consider molecular ly 
repelling pores where < a < ao (Ea > 0, Ep = 0); 
here, p « eKp(—Ep/kBT), independent of the pore en- 
ergy level. As the pore is made attracting, p will acquire 
ex.p{Ei3/kBT) behavior and can be defined as 



P = Po 



—9{a — ao) + d{ao — a) 

"0 



(9) 



where 9{x > 1) ~ 1 is the Heaviside function indicat- 
ing when the pore first becomes molecularly attracting. 
Using (|^), J{N^a — > 00) (the current in the limit of in- 
finitely attracting pores) becomes 



J°°(7V) 



ao(V-l)(l- A)+po(2-A)' 



(10) 



which can approach A/(2 — A) > J*. For ^ a*, the 
maximum remains at 0* . However, when ap a* , a, 
p « exp^Efj/kBT) and the maximum in flux as a func- 
tion of a is preempted by a current which approaches 
(|^). The condition for J°°{N) < J*{N) (a remaining 
maximum in J as pore well depth is increased) is deter- 
mined by 



r,*2 



ao > a 



1 



(11) 



Fig. ||(b) compares the behavior of 10* x J(5) using 
P = Po — 0.1 (solid curve) with that of 10* x J(5) using 
Eq. d) (broken lines). For A = 0.02, the maximum at 
a* = a/15/14 is destroyed when ( [ll| ) is satisfied, ao ^ 
1.562. Estimating ao from (||) and A = 0.02, ao < a* 
for gases, but can be 0(1) for liquids. Curve PqPiP^^P^ 
retains the maximum J* since ao — 10°'^^ > 1.562, while 
PaPiP2 corresponds to ao = 10^"'^^ < 1.562 which gives 
a monotonic J as a ^ 00. A high current may occur 
at a — > 00 despite the high pore occupancy due to an 
accompanying exponential increase in p. As expected, 
these high values are more difficult to achieve as ao (as 
well as N) increases, because the onset of exponentially 
increasing p is delayed. 

Now consider the nonlinear effects of large A. Accord- 
ing to (|^), values of a that give a maximal J* depend 
strongly on A; thus, A ^ 1 can yield large a* » ao 
where the maximum in flux is destroyed. At the max- 
imum value A = 1 (corresponding to approximately to 
pure solute or vacuum in (i?)), the maximum flux occurs 
when p ^ N and approaches J(A = 1) ~ a/(a 4- 2). 
The smallest flux occurs when p <^ N and a <C 1 as 
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expected. The nonhnearity of J(A) is important only 
when a ^ 1, as shown in Figure ||, corresponding to a 
pore interior with high particle occupation, when particle 
exclusion nonidealities are most pronounced. 

p=LO, N=10 
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FIG. 3. Nonlinear behavior of 10^ x J(10) for p = 1.0 as 
a function of a. Note the competition between nonlinearities 
in a and A, particularly at large A. 

To experimentally test possible nonlinearities densities, 
pressures, and temperatures may need to be substantially 
varied, and is readily performed only in zeolite systems. 
Adiabatic ultrasonic driving with frequency oj <^ J^^ of 
the fluid in (R) (not affecting a which depends only on 
Ul) will enhance transport: Upon setting A{t) = Aq + 
Ai cos Lot and averaging over a period, 

^ J(A^,Ai=0) 
J(iV) = J{N, Ai = 0) + ^ ' ' X 



p « (a + l)iV 



(12) 



0(Ai)3 



for small Ai/Sq. In the p:s> {a + 1)N case, the tempera- 
ture dependence of the flux enhancement is —iEa/ksT 
or —Ea/ksT, while the second case is temperature inde- 
pendent. 

Effects of microscopic internal pore structure/disorder 
on microflow can also be modelled. Consider N* defects 
in hopping rates p^(i) = ^^(Oi between sections i{k) and 
i{k) + l. When pl{t) is statistically independent of n^t), 
the mean-field result 



J 



p{aP — 7(5) 



{a + j){l3 + S){N -l + n)+p{a + (3 + j + S)' 

(13) 



where il = ^k=i(p/P*k ~ 1)' exact. Biological and 
zeolitic channels contain internal pore structure expos- 
ing random binding sites Ep(k). An example of a one- 
defect pore is gramicidin A channels across a biomem- 
brane. These are composed of two barrel structures, 



one in each lipid layer, joined by attractive molecular 
interactions near the bilayer midplane forming a defect 
bisecting the entire channel. Osmosis experiments on 
Gramicidin A/bilayer liposomes reveal rich temperature 
dependences which are interpreted as lipid phase tran- 
sitions inducing changes in how the Gramicidin A bar- 
rels are joined |l^. Molecular permeation through pure 
lipid bilayers is also modelled by (|l^). Biolipid head- 
groups offer activation barriers such as those depicted 
in Fig 1^. Within the hydrophobic lipid tails, particle 
motions are predominantly perpendicular to the mem- 
brane along the aliphatic chains. Thus, bilayer mem- 
branes physically resemble many close-packed channels 
with defects at the bilayer midplane, where the ends 
of the hydrocarbon chains approach each other, and 
at stiff unsaturated bonds along the fatty acid chains. 
Various zeolites are comprised of interconnected cages 
and joints which are also modelled using (^3|), although 
actual flow measurements are ensemble averages over 
macroscopic membrane regions containing many pores: 
{J) = 1:n,N',p' .f{N,N*,pl)JiN,N*,pl), where d is the 
membrane thickness, Nmin — d/£, and f{N,N*,p'^) is 
the distribution of channels with arc-length L = N/i, 
number of defects N*, and defect strengths p^. We men- 
tion that the steady state result (|lj) is also an exact 
solution to the Heisenberg spin chain hamiltonian with 
boundary conditions determined by {/i} and quenched 
random energies ]l^ . 

Since the particle density at section i obeys = 
Qiji+i ^ + Ti_i) which on scales 3> ^ is a diffusion 
equation with qi^ defining a cooperative diffusion co- 
efficient, osmosis occurs strictly by mass diffusion, dif- 
ferent from tracer diffusion. Gonditions required for a 
maximum in J(a) depend on the pre-exponential factor 
ao for realistically defined {/i}. In LTE, relationships 
among the various kinetic parameters and various fluids 
can be determined by equilibrium measurements such as 
solvent-solute volumes and heats of mixing. The different 
Arrhenius temperature regimes delineated may also be 
simple way of systematically probing solvent-pore inter- 
actions. Experimentally the nonlinearities in A may be 
probed by low frequency acoustic driving. Electrostric- 
tion and mechanical deformation of the pores may also 
control a and p via Vp Q . Effects of unstirred layers near 
the pore mouths can be easily treated with macroscopic 
convection-diffusion equations, but are experimentally |^ 
and theoretically fl^ found to be unimportant in aque- 
ous osmosis. When important, unstirred layers yield an 
implicit equation J oc A (J) ||l^. The simple model pre- 
sented along with the numerous applicable physical sys- 
tems and proposed experimental tests complements MD 
simulations and will subsequently lead to a better under- 
standing of more complex systems, including multispecies 
transport and chemical reactions in pores. 
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